Continuous Avalanche Segregation of Granular Mixtures in Thin Rotating Drums 
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We study segregation of granular mixtures in the continuous avalanche regime (for frequencies 
above « 1 rpm) in thin rotating drums using a continuum theory for surface flows of grains. The 
theory predicts profiles in agreement with experiments only when we consider a flux dependent 
velocity of flowing grains. We find the segregation of species of different size and surface properties, 
with the smallest and roughest grains being found preferentially at the center of the drum. For a 
wide difference between the species we find a complete segregation in agreement with experiments. 
In addition, we predict a transition to a smooth segregation regime — with an power-law decay of the 
concentrations as a function of radial coordinate — as the size ratio between the grains is decreased 
towards one. 



Mixtures of grains differing in size, shape or density 
segregate in vibrated containers Q, in two and three- 
dimensional rotating drums and by pouring the 
mixture into a heap For instance, axial band segre- 
gation appears in long rotating cylinders filled with mix- 
tures of grains |^-^| , while radial segregation of granular 
mixtures is observed in thin (w 10 grain size) rotating 
drums: large grains are usually being found near the 
outer region o_f the drum, and small grains near the center 
of the drum 



Theoretically, the segregation of mixtures in rotating 
drums has been studied by Zik et al. || using a contin- 
uum approach based on mass conservation to treat band 
formation in 3-D drums. Recently, Khakhar et al. Q pre- 
dicted the concentration profiles in 2-D rotating drums 
filed with mixtures composed of equal-size grains of dif- 
ferent density. However, there is no general analytical 
solution for segregation of grains in 2-D rotating drums 
in the case of mixtures of arbitrary size, and shape or 
friction properties. In this Letter, we study analytically 
this case focusing on the continuous avalanching regime, 
where a steady state is observed. Based on a continuum 
theory of surface flows of grains, we obtain predictions 
for the profile of the rolling species which agree quantita- 
tively with experiments when we take into account 
the velocity profile of the flowing grains. 

We find that the small and rough grains occupy the 
center of the drum while the large and rounded grains 
segregate to the bottom of the drum. According to the 
size ratio and the different angles of repose of the species, 
there are two distinct regimes. For large size ratios our 
solution shows almost complete segregation of the mix- 
ture with a small region of mixing near half of the drum 
radius. The surface concentrations of grains behave ex- 
ponentially on radial coordinate and segregation is char- 
acterized by a typical length scale inversely proportional 
to the size ratio between the grains. These findings agree 
with the experiments of Clement et al. Q . On the other 
side, when the species differ slightly, we predict a slow 
power-law decay of the concentrations of grains at the 
surface, with no characteristic length scale associated to 
segregation. 



In the continuous avalanching regime, grains on the 
free surface flow down steadily due to gravity forming a 
film of rolling species, while the bulk of the material ro- 
tates as a solid body. Thus, it is plausible to describe the 
system only by its surface properties ||l0f-|l3||. 

We consider two local coarse-grained equivalent thick- 
nesses of the species in the rolling phase Ri(x,t), with 
i = 1, 2 [i.e. the total thickness of the rolling phase mul- 
tiplied by the local volume fraction of the i grains in 
the rolling phase at position x on the surface of the bulk, 
Fig. 0(a)]. The total thickness of the rolling phase is then 
R(x,t) = Ri(x,t) + R 2 (x,t). The static phase rotating 
as a bulk is described by the local angle of the surface 
height, 9{x,t) = —dh{x,t)/dx [h(x,t) is the height], and 
the surface fraction of static grains, 4>i{x,t), of type i at 
the surface of the bulk. The equations of motion for the 
rolling species are |L2| 



dRi(x,t) 
Ft 



v(x, t 



dRi 
dx 



(0.1a) 



and the surface concentrations of static grains are given 
by 



,.dh 



wx) = -Yi 



(0.1b) 



and 4>\+ <j>2 = 1- Here w is the rotation frequency of the 
drum. The downhill convection velocity of the rolling 
species v(x,t) is usually taken as a constant P,[lO| [l3| . 
However, experiments on granular flows in drums, in- 
clined chutes and sandpiles (|jl^,|l^] indicates that the 
velocity varies linearly with the vertical position, y, of 
the particle in the flowing layer, i.e. v(y) = 2ay, where 
a is the shear rate. In our continuum approach we take 
a height-averaged velocity dependence of the form |l5| : 
v{x,t) — 1/R(x,t) J Q v(x,y,t)dy — a R(x,t). We will 
show that the flux dependence of the velocity plays a 
crucial role in explaining the experimental results. 

The interaction term Ti takes into account the conver- 
sion of static grains into rolling grains and vice versa, and 
is @, [i= 1,2, (JjLi)] 



Yi 



R l ~b l {e)R l +x l {B)4> l R J 



(0.2) 
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The interaction term is a function of the collision func- 
tions contributing to the rate processes. Amplification, 
ai(9): an i static grain is converted into a rolling grain 
due to a collision by a rolling grain i. Capture, bi{9): 
an i rolling grain is converted into a static grain. Cross- 
amplification, Xi (0) : the amplification of a j static grain 
due to a collision by an i rolling grain. 

We consider the geometry of a two-dimensional rotat- 
ing drum of radius L [-L < x < L, Fig. [I] (a)]. We study 
the continuous avalanche regime where steady state pro- 
files are observed in the drum, and we seek a solution 
dh(x)/dt — 0, dRi(x)/dt — 0, with boundary conditions 
Ri{x = -L) = Ri(x = L) = 0. 

The profile of the total thickness of the rolling phase 
R{x) and the thickness of the rolling layer at the center 
of the drum A = R(x = 0) are obtained from (0.1): 



R{x) 



(L 2 



A = R(x = 0) 



r 2Nl/2 



L 



(0.3a) 



(0.3b) 



The square root dependence in (0.3a) differs from the 
parabolic profile found in (ll]]|] . This latter result [fllpl 
is obtained by assuming v =const. Solving Eqs. (0.1) 
for this case gives: R(x) = w/(2v)(L 2 — x 2 ), and A = 
wL 2 /(2v). In Fig. |2|(a) we plot X(w) and compare with 
the digital photography experimental data of Khakhar et 
al. H for spherical sugar balls of 1.8 mm size, and the 
MM experiments of Nakagawa et al. . We clearly see 
that X(w) is con sistent with a square root dependence 
on w as in ( |0.3b ) and not with a linear dependence on w 
obtained by assuming v =c onst fl6||. In Fig. 0(b) we plot 



the theoretical prediction (3.3a) for R(x), as well as the 
prediction obtained when v =const, and compare with 
the experimental results of S. Since we use a — 25/s 
obtained from Fig. §(a), and w — 2nf with / = 3 rpm 
from the experiments, the comparisons in Fig. ^](b) are 
made without tuning of parameters. We again see an 
excellent agreement between experiment and theory only 
when the flux-dependent velocity is taken into account, 
while the prediction based on the assumption of a con- 
stant velocity profile gives lower values of R(x) ]l7| . This 
result is in agreement with the numerical analysis of j^] 
for simple shear flows. 

To describe the system composed of a binary mixture 
of grains, we propose below a set of collision functions 
valid in the limiting cases of a wide and small difference 
between the species. We first treat the case where the 
difference in size and angles of repose is not too large 
(we will quantify this later). In this case the region of in- 
terest is concentrated in a small region around the angles 
of repose, a fact that allows us to linearize the collision 
functions in the vicinity of the angles of repose. We pro- 
pose [Fig. (b)] 



cn(9)=Xi(fi) = C + 7 [9{x) 
h(9) = C~ 1 [9(x) 



(0.4) 



Here, 7 ss 23/s |15) and C are collision rates. In the 
continuous avalanche regime we have w > 7<5 2 /2, where 
6 is the difference between the maximum angle of stabil- 
ity and the pure angle of repose of the species Jll[]. We 
note that the collision functions (0.4) differ from the ones 
proposed in jl^] in the crucial fact that the generalized 
angle of repose, 9i(<j>j), of a rolling grain type i depends 
on the composition of the surface <pj 101 : 



9i(</> 2 ) = m4> 2 + 
h{4>2) = W>2 + 021 



-m<j>i + 022- 



(0.5) 



The limiting cases are 0jj = 9i(cj>j = 1), m = $12 — 0n = 
022 ^ 021, and the difference rp — 0i(02) — 02(^2) is inde- 
pendent of the concentration 4>2, then tp = 9\\ — 02i = 
012 — 022- 

The angle ip is proportional to the degree of differ- 
ence in size and shape between the species, and there- 
fore determines the degree of segregation; the larger t/j, 
the stronger the segregation. If the species 1 are the 
smallest then 0i(02) > 02(^2) for any composition of the 
surface ^2, and also we have 0i2 > 021- For mixtures 
of grains with different shapes or friction coefficients we 
have 0ii 7^ 022, and 0i2 = 02i if the species have the same 
size. If 0ii is the repose angle of the pure round species 
and 022 is the repose angle of the pure rough species, then 
0n < 022- If the species have the same shape or friction 
coefficients then 0n = 022- When the size ratio is close 
to one (d 2 /di < 1.5 jUP, rfj size of grain i), the angle ip 
is s mal l and then we linearize the collision functions as 
in (0). 



Using (0.1b) and (0.4) we obtain the profile of the sur- 



face as a function of the rolling species 



9{x) - 9^) 



-4>R2 
R 



(0.6) 



where we use that di/ipx <C 1. The equation for the 
rolling species of type 1 is obtained from Eqs. (0.1) and 
& 



dRi(x) dR{x) 



dx 



dx 



■07-R2 
CR 



Ri 
R 2 ' 



(0.7) 



By setting u(x) = R\[x)/R{x) we obtain u'/[u(l — u)] — 
—(; x/[L 2 — x 2 ], where £ = jtp/C, and we arrive to a 
power-law form 



Rx(x)/R(x) = 



1 



1 + A(L 2 



(0.8) 



where the power-law exponent ^ ~ ip w 0.2 |T^] de- 
pends on the degree of difference between the species, 
and A — ((f> 2 /$i)L 2 ^ is an integration constant obtained 

from (1/2L) J^ L (f>i(x)dx = <I>i, where $i is the fraction 
of species i at the surface (usually <I>i = 1/2). The surface 
concentration profile shows a slow power-law decay 



2 



0i(a;)=u[l + £(l-u) 



(0.9) 



Figure || shows the solutions obtained for this case (small 
ip). The species 1 — the smallest grains or the rough- 
est grains — are found preferentially at the center of the 
drum [0i (x = 0) > <p2 (x = 0)] . The degree of segregation 
is small; the surface concentrations behave very smoothly 
as a function of the position in the drum. 

According to (3.9), the exponent £, controls the inten- 
sity of segregation. Since £ ~ ip oc d^/dx, then as we 
increase the size ratio, ip increases, and the linear ap- 
proximation of the collision functions (0.4) breaks down. 
In this case, strong segregation effects act in the sys- 
tem. Segregation happens in the flowing layer as the 
large grains are observed to rise to the top of the rolling 
phase, while the small grains sink through the gaps left 
by the large grains; an effect known as kinematic siev- 
ing, free-surface segregation or percolation [|l8|,[l5| . Then 
the small grains are captured first near the center of the 
drum, since they interact with the surface before than 
the large grains. Thus, for large if) (in terms of size ratios 
(fo/di > 1.5 H]) the capture function of the large grains 
& 2 vanish near the angle of repose of the small type 1 
species: the small rolling grains screens the interaction 
of the large rolling grains with the surface. We simulate 
this effect by considering the following non-linear colli- 
sion functions (]]J (valid for large ip): 



h{9) ee 7 n[0i(^) - 0(a;)], 



(0.10) 



where TL[z] = 0, if z < 0, and U[z] = z, if z > [Fig. | 
(c)]- 

Next we find the solutions for x < using (3.10) ||l9|. 



We consider that there are R® and i?2 rolling species 
of type 1 and 2 res pectively at the center of the drum 
(Rl + R% = vWai, and E%/{R\ + JJ§) ~ $ 2 ). Since 
the collision functions are defined according to the value 
of 6(x) in comparison to 0,(02), we divide the calculations 
in two regions: inner region, where 62 (02) < 9 < 0i(02), 
and outer region, where 9 < 9 2 ((p2) < #i(02)- 

Inner Region. If 62 < 9 < 9%, we obtain from ( 3. lb| ) 
and (0.10): (pi(x) — l,cj)2(x) = 0. Then we obtain the 



profile of the rolling species using (0.1) 



Ri(x) 



™(L 2 
a 



x 2 ) 1/2 -R° 2 , R 2 (x)=R%. (0.11) 



The profile of the surface is: 
This solution is valid when 
6*21, or for > x > 



9(x) 
9(x) 



> 



= wxRi{x)/^. 
6 2 {<t>i = 1) = 
[&2y/ wa /l' l P ~ 



Vl-*|+ wa/(jip) 2 ]L/[l + wa/(^) 2 }. 

The profile of the small rolling grains decays with the 
distance fro m the center as in the case of the total rolling 
species Eq. (0.3a), since the small grains are the only one 
effectively interacting with the surface due to percolation, 
so that it is a single species problem. Next, we find the 
solution of the problem at the lower end of the drum. 



Outer R egion . If 8 < 9 2 < 9\ {—L < x < x m ), using 
(0.1) and ( 0.10| ) we obtain the equations for the rolling 
species 



dRi{x) 
dx 



(iVX — r )1pR2) 



Ri 



aR{xf 



(0.12) 



To solve this equation we set u = R\{x)/R{x), and ob- 
tain u' /[(l — u)u] = jip/[aR(x)] and the solution is 



R!(x)/R(x) 



1 



1 + Aexp[— (x — x m )/a 



(0.13) 



where the characteristic length of segregation is x s = 
■\/waL / (jip) ~ di/tp, A — $ 2 /®i is an integration con- 
stant obtained by considering the continuity at x — x m , 
and we have replaced exp{sin _1 [(a; — x m )/L]} ~ exp[(x — 
x m )/L]. The concentration profiles is 



<pt(x) = exp[(x - x m )/x s 



(0.14) 



Figure || shows the theoretical profiles found for this 
case (large ip). At the center of the drum we have to- 
tal segregation of the small grains [4>i(x = 0) = 1, 
02 (x = 0) = 0]. Then the concentration of small grains 
decays exponentially — with a region of mixing of char- 
acteristic length x s . We also notice that in this case 
the profiles depend o n w [x s ~ y/w] while in the case of 
small ip, the profiles (0.9) are independent of the rotation 
speed. Typical values of the phenomenological constants 
of the theory have been obtained experimentally in [ p"5"| : 
7 w 23/s, a w 25/s [see also Fig. ||(b)], an d V> ~ 0.2. 
Then we obtain x s ~ 1.6 cm. These results are in quan- 
titative agreement with experiments of Clement et al. || , 
who found an exponential decay in the concentration of 
particles with a characteristic length of sa 1.7 cm for size 
ratio d 2 /di = 1.5. We also notice that x s decreases as 
d 2 /di increases, since ip oc e^/di; a relationship found in 
the experiments of @, too. 

In sum, we show that the predictions based on the com- 
mon assumption of a constant velocity profile of flowing 
grains are in error. The profiles predicted by the the- 
ory agree with experiments only when a height-averaged 
velocity profile is taken into account. If the difference 
between the species is wide, we find a sharp exponential 
segregation in agreement with experiments. It would be 
interesting to test experimentally the transition to the 
length scale- free segregation regime when d 2 /d\ is de- 
creased below ~ 1.5. We believe that the present anal- 
ysis has the potential to unify the mechanisms underly- 
ing segregation in different geometries such as silos and 
rotating drums. Our results are also valid for granular 
mixtures poured in 2-D silos: the different segregation 
regimes have been observed in [n2|Jl3|, and experimen- 
tally by Grasselli and Herrmann B. 
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College de France for hospitality in initial stages of this 
work. 
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FIG. |l](a) The thin rotating drum, (b) Linear approx- 
imation (solid lines) to the collision functions (dashed 
lines) when the difference between the species is small, 
(c) Non-linear approximation (solid lines) to the colli- 
sion functions when the species differ appreciably. The 
collision functions are expected to be continuous (dashed 
lines), but we approximate them as discussed in the text. 

FIG. |(a) Theoretical prediction ( |0.3b| ) for X(w) and 
comparison with experiments of 00]- (b) R(x) profiles 
from theory and experiment || . We also plot the predic- 
tion obtained with the assumption v =const. 



FIG. H Theoretical profiles <fii(x) for small (dashed 
lines) and wide (solid lines) differences between the 
species (small and large ip, respectively). 
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